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Abstract. A formulation of Continuum Mechanics within the context 
of General Relativity is presented that allows for the incorporation of 
certain types of anelastic material behaviour, such as viscoelasticity and 
plasticity. The approach is based on the concept of a four-dimensional 
body-time complex structured as a principal bundle. The anelastic re- 
sponse is regarded as the result of a continuous distribution of inhomo- 
geneities, whose evolution is dictated by a suggested relativistic version 
of the Eshelby tensor. The role played by various groups is emphasized 
throughout the presentation and illustrated by means of the example of 
an anelastic fluid. 



1. Introduction 



There are many situations in physics where a knowledge of the thermo-mechanical 
ly-j ' bulk properties of matter is essential for an understanding of a wide variety of nat- 

ural phenomena. Among these properties the response of matter to internal and 
' external forces is often a dominant feature that determines its static or dynamic 

q-i behaviour. The laws of classical (i.e., non-relativistic) Continuum Mechanics have 

j 1 ^ ' been developed to address such phenomena and constitutive theory is now a ma- 

(^j), ture branch of material science. The formulation of traditional non-relativistic 

Continuum Mechanics is firmly rooted in the framework of Newtonian dynamics 
and time dependent transformations that preserve the Euclidean structure of a 
' three dimensional space are deeply ingrained in its insistence on the principle of 

material objectivity 22 and compatibility with thermodynamics. Furthermore, the 
evolution and past history of material deformation depend on a universal New- 
tonian time and are therefore incompatible with notions of observer dependence 
in a relativistic spacetime. Whilst such discrepancies may be irrelevant for some 
problems, they should not in principle be ignored. They cannot be ignored for 
problems in which relative material speeds comparable with that of light are rele- 
vant or for problems involving non-Newtonian accelerations or large gravitational 
tidal stresses. Such situations can arise in numerous problems in astrophysics and 
space science. Modelling the behaviour of rapidly rotating neutron stars offers a 
challenging avenue for exploring new material phenomena as do situations involving 
the interaction of matter with intense electromagnetic or gravitational radiation. 
The constitutive properties of materials in such situations are far from clear since 
they depend on particular reformulations of the laws of non-relativistic Continuum 
Mechanics that must be compatible with the local Lorentz structure of spacetime 
needed to maintain the causal propagation of physical effects within all media. 

A number of relativistic reformulations exist. Based on pioneering efforts by 
Rayner |2T] . Carter and Quintana 0, Maugin Kijowski \T^\, and others, 



2 



MARCELO EPSTEIN, DAVID A BURTON, AND ROBIN TUCKER 



most reformulations model the spacetime histories of a three-dimensional mate- 
rial medium by a congruence of time-like worldlines in spacetime M. and identify 
a material body B with a class of three-dimensional spacelike surfaces transverse 
to such a congruence. It is assumed that equations can be prescribed to determine 
a map / : M. — > B such that in any coordinates X A for some domain of B the 
components f A of / satisfy d f 1 A d f 2 A d f 3 ^ 0. Using the spacetime metric 
g, the 3— form d f 1 A d f 2 A d f 3 is mapped to a future pointing time-like normal- 
ized 4— velocity field, the integral curves of which model the world lines of material 
elements of B in M.. The equations for / can be derived from a variational ac- 
tion principle or (for dissipative phenomena) by setting to zero the divergence of a 
suitable stress-energy tensor on M . The applications of this formalism are largely 
dictated by how the action or stress-energy tensor are constructed in terms of / 
and the spacetime metric g. Since, in general, spacetime has no particular Killing 
isometries the role of the Euclidean structure in non-relativistic Continuum Me- 
chanics is at best a guiding principle in the modelling process and a new principle 
of objectivity must be found within covariances of the spacetime formulation. 

The dissipative behaviour of many materials can be modelled within the context 
of classical (i.e., non-relativistic) Continuum Mechanics by means of an elastic 
prototype, namely, a fixed elastic material "point" . The most striking example of 
this situation is provided by the theory of ideal elasto-plasticity, whereby the process 
of plastification does not alter the elastic constants of the material, as evidenced 
by subjecting the material to an unloading process. What does change during 
plastification is, therefore, the way in which the clastic prototype is implanted 
into the body as time goes on. A similar interpretation can be given to processes 
of biological bulk growth, such as those observed in bone and muscle. Here, the 
anelasticity reflects the fact that, as time goes on, more (or less) material of the 
same type is squashed into the body, the properties of the new material being 
identical to those of the existing one. As in the case of plasticity, the result is, 
in general, the development of deformations and of residual stresses. We may 
say that in these phenomena there are two different types of kinematics at play: 
the ordinary kinematics of embeddings of the body manifold into space and the 
material kinematics resulting from the time variation of the implants. These two 
kincmatical paradigms are quite independent from each other. Correspondingly, 
the forces behind them are of diverse natures. As much as ordinary (Cauchy or 
Piola) stresses can be seen as the causes of the evolution of the spatial kinematics, 
the forces reponsible for the material evolution are the so-called configurational 
stresses, pioneered by Eshelby as far back as 1951 in his classical paper ^5] dealing 
with the force acting on an elastic singularity. Assuming we start from an originally 
perfectly homogeneous stress-free body and we subject it to some loading process, 
if the material is of the type just described, the different points may undergo a 
material evolution (plastification, say), each point to a different extent than the 
others in its neighbourhood. The result (after unloading, for example) is that, 
although the body is still uniform in the sense that it is made of the same material 
at all of its points, it can no longer be considered homogeneous: it has an irreversible 
continuous distribution of inhomogeneities in its midst. We may say that this type 
of inhomogeneity and the previously described type of anelasticity are the material 
and spatial counterparts of each other. 
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In any relativistic framework it becomes necessary to accommodate material 
features that have a natural interpretation in a non-relativistic limit in terms of 
classical hysteretic dissipation, plasticity material inhomogeneity, aging and more 
general anelastic properties that involve memory effects and rate dependent re- 
sponses to the environment. In this paper we begin an exploration of these issues 
in terms of a new formulation of relativistic Continuum Mechanics based on the con- 
cept of a four-dimensional body-time complex Af rather than the three-dimensional 
body manifold B. The essential idea in our relativistic formulation is to use a rank 
4 map K : Af — ► A4 to describe the physics of a deformable body. Embeddings of 
lower rank have been used to define relativistic point particles, strings and mem- 
branes [HH, |23, [El, ESI, d and for higher dimensional Af and M (and 
reparameterisation invariant actions), which form the cornerstone of current efforts 
to model the basic interactions. It is our intention to base a general relativistic 
theory of constitutive properties of materials on the map k that offers a simpler 
approach for a relativistic description of anelasticity. 

2. The material body-time manifold 

In classical Continuum Mechanics, the notion of a material body B as a three- 
dimensional differentiable manifold can be said to emerge from the following consid- 
erations. All possible spatial manifestations (or configurations) of an identifiable 
entity happen to be connected open sets in R 3 with the property: there exists 
a unique and smooth correspondence between the points of each and every pair 
of configurations. We call the corresponding equivalence classes material points. 
Since the correspondence has been assumed to be smooth with respect to spatial 
positions, it follows that each configuration can be seen as an atlas of an under- 
lying differentiable manifold B, called the material body. It is common to assume 
that there exists an atlas consisting of just one chart. The configurations can now 
be regarded as embeddings of the material body in M 3 . This fastidious heuristic 
"derivation" of the notion of a classical material body will serve the purpose of 
justifying the following relativistic counterpart. 

In General Relativity the object under consideration manifests itself always as 
a world tube T consisting of a congruence of world lines in a relativistic four- 
dimensional space-time manifold AA. The assertion that this is a congruence of 
world lines and not just a world tube is a reflection of the fact that, just as in 
the classical counterpart, we postulate that particles are distinguishable from each 
other. Moreover, on each of these lines there exists a unique-up-to-translation 
parametrization such that: 

(2.1) g(u,u) = -l, 

where g is the metric of AA and u is the future-pointing timelike unit vector field 
tangent to the world line. In each of these tubes, therefore, there exists a physically 
meaningful (local) action of the group of translations of K and, for the sake of 
brevity, we denote this group also by M. The manifold of integral curves of u is the 
quotient B = T/K and we are naturally led to the following definition: 

Definition 2.1. A material body-time (or body-time complex) is a principal fibre 
bundle (Af, B, ir, K, R). The total manifold Af alone will be referred to as the ma- 
terial body-time manifold. The base manifold B will be assumed to be an ordinary 
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material body, namely, a trivial three-dimensional manifold. The projection map: 

(2.2) tt:N — >B, 

assigns to each body-time point the corresponding body point. The principal nature 
of the bundle reflects itself in that the typical fibre M and the structural group 
coincide. 

A material tangent vector V £ Tfif is said to be vertical if 

(2.3) tt*(V)=0. 

We endow the fibre bundle (TV, B, ir, E, R) with the fibre metric g R = —dS (g> dS 
where S is a fibre coordinate and define on A/" a vertical vector field U = ^ 
satisfying 

(2.4) g R (U,U) = -l. 

3. Histories 

A history is the relativistic counterpart of a classical configuration. As discussed 
above, therefore, a material history can be identified with a congruence of world 
lines. More formally, we state the following: 

Definition 3.1. A history is an embedding 



(3.1) k:JV — >M, 
such that: 

(3.2) g(K*(U), K *(U)) = -l. 

Since U = ^ Equation (|3.2I) identifies S with the proper time of each integral 
curve of k*(U). The tangent map 

(3.3) : TAf — >TM, 



appearing in Equation lj3.2|l is called the history gradient of k. At each point n G Af 
it is represented by a tensor 1 F = k» : T n Af — > T K ( n )M. 

Remark 3.1. For simplicity, we are assuming that a history is an embedding of 
the whole of Af. In reality, one may have to consider histories that consist of 
embeddings of some time- wise portion of Af. 

We will often refer to Equation (|3.2I) as the time consistency condition of our 
embeddings. It represents a fundamental constraint in the dynamical description, 
as we shall soon see. 

Remark 3.2. It would also have been possible not to require the satisfaction 
of the time-consistency condition. In other words, one can derive an alternative 
theory whereby the structural group of Af is the group of all (orientation-preserving) 
diffeomorphisms of the real line. It is not entirely clear at this point whether both 
formulations are equivalent up to a re-parametrization of solutions. This may very 
well be the case under some constitutive restrictions (some of which are discussed 
in later sections), but it appears that in the most general situation the consistency 
condition (and its associated Lagrange multiplier) play a fundamental role in the 
theory to be presented. 



Since the vector spaces connected by this tensor are associated with different manifolds it is 
often referred to as a two-point tensor 
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It may prove convenient to exhibit the local coordinate versions of some of the 
basic quantities and equations introduced so far. All coordinate systems in AT will 
be assumed to be consistent with a trivialization of the fibre bundle. The four 
coordinates, therefore, will be split into three body-coordinates X A , A — 1,2,3, 
and the proper time-coordinate S. For compactness, we will also often adopt the 
notation X A , A — 1,...,4, with X A = S. In a spacetime domain, on the other 
hand, we will denote the coordinates by lower-case letters: x a , a = 1, ...,4. Thus, 
a history n is given in coordinates in terms of four smooth functions: 

(3.4) x a = K a (X A ) = K a (X A ,S), 

and the corresponding history gradient is represented by the field of matrices: 

(3.5) n = *£. 

This matrix field is of maximal rank (4) , since we have assumed the histories to be 
embeddings. We will also assume it has a positive determinant. The image k*U of 
the canonical unit vector field U = ^§ provides us with the four-velocity field u on 
M, given in components by: 

(3.6) u a = F%U A . 

With M. time-orientable u will be declared future-pointing for S > 0. Since the only 
non- vanishing component of U is U 4 — 1, condition l|3.2|l becomes in components 

(3.7) gabF^Fl = -1. 

4. Elasticity 

In classical Continuum Mechanics, a first-grade (or simple) material point with- 
out memory effects is characterized by constitutive laws (for the stress, the heat 
flux, the internal energy and the entropy) whose independent variables are the 
present local values of the deformation gradient, the temperature and the temper- 
ature gradient. A material is said to be elastic if the stress tensor depends on the 
deformation gradient alone. Thus, in the classical setting, it is possible to remain 
within the realm of a purely mechanical theory, whose field equations are the con- 
servation of mass and the balance of linear (and angular) momentum and whose 
constitutive response is just an elastic law for the stress tensor. In a similar vein, 
one may avoid a complete and explicit thermodynamic treatment even in a theory 
of viscous flow, by allowing time-derivatives of the deformation gradient to enter 
the picture, but ignoring the temperature and its gradient. If an elastic material 
is subjected to an internal constraint, such as incompressibility or inextensibility 
along a material direction or any other restriction imposed on its deformation gra- 
dient, the stress is not completely determined by the deformation gradient, since 
the forces (stresses) necessary to maintain the constraint are not related to the 
deformation. For example, if a material die is incxtcnsiblc in the direction of one of 
its sides, one may apply an arbitrary surface traction in the inextensible direction 
without producing any deformation! The total stress is then assumed to be resolv- 
able into the sum of two parts: one determined by the deformation (according to 
some elastic law) and another one given by a Lagrange multiplier applied to an 
expression dictated by the constraint equation. The standard way to obtain this 
expression (see [221, Section 30) is based on the assumption that the second part 



6 



MARCELO EPSTEIN, DAVID A BURTON, AND ROBIN TUCKER 



does no work in any motion satisfying the constraint. In other words, the constraint 
forces do not expend mechanical power on the admissible virtual velocities. In the 
case of holonomic constraints, this criterion leads to the gradient of the constraint 
equation. More precisely, let <j)(F) = be the given constraint equation, where 
F is now the classical deformation gradient. Then the term to be added to the 
determinate part of the stress is of the form: 

(4.D A§F- 

A being the Lagrange multiplier. 

In General Relativity the thermomechanical response of a material is assumed 
to be encoded in the symmetric stress-energy-momentum tensor t. In trying to 
emulate the classical theory of elastic simple materials, one can assume that the 
value of t at a point m of M. occupied by a point n G J\f depends exclusively on 
n, on the local value (at m) of the metric g and on the local value (at n) of the 
history gradient F. This a-priori exclusion of thermal effects may need some further 
discussion. Be that as it may, it is important to realize that in our kinematical 
picture all materials are to be considered as internally constrained. The reason 
for this somewhat surprising assertion lies in our time-consistency condition l|3.2[) . 
which is a (universal) restriction on all possible history gradients. Following the 
paradigm of classical Continuum Mechanics as given in Equation Q4.1|) . we obtain 
the following general form of the stress-energy-momentum tensor of a relativistic 
elastic material: 

(4.2) t = f (n, g, F) + A(n) u <g> u, 

where, as before, u = FU, and where f is an arbitrary symmetric tensor-valued 
function of the arguments shown. The particular form of the term affected by the 
Lagrange multiplier A in (|4.2[1 arises when applying the prescription (|4.1|l to the 
constraint equation i|3.7[l . namely with: 

(4.3) cj> = g ab F%F b B U A U B + 1. 

Note that Equation l|3.7fl is imposed (i.e. is set to zero) after calculating Q4.1fl , 

An elastic constitutive law is said to be hyperelastic if f derives from a scalar 
potential W = W(n, g, F) according to the formula: 

(4.4) f = Wg -i + ^pT 
or, in components: 

, , ,dW u 

(4.5) f ab = Wg ab +g *dJ_ F b A . 

The quantity W corresponds to the classical strain energy per unit spatial volume. 
In general relativity, therefore, we interpret it as the elastic potential per unit space- 
time volume as induced by the metric, namely by the 4-form in A4 with component 
y/— det[g a b] in a coordinate basis. An alternative approach to determine the gen- 
eral form of the hyperelastic stress-energy-momentum tensor is via a constrained 
variational principle, as shown in the Appendix. 

Remark 4.1. In classical Continuum Mechanics it is customary to introduce a so- 
called reference configuration and to induce the Cartesian volume therein onto the 
body manifold. In that case, it makes sense to define a strain energy density per unit 
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volume in the body (but this volume depends of the reference configuration chosen 
through the determinant of the gradient of the change of reference) . Alternatively, a 
volume form is introduced in the body representing the mass density, in which case 
one can introduce a strain-energy per unit mass. By defining various pulled-back 
stress tensors or tensor densities, the expression (|4.4ll is simplified to the extent 
that it loses its "spherical" part, retaining only the term with the partial derivative 
of the potential (see [22] , Section 82). Similar arguments can be used in general 
relativity, the most widespread being that of the baryon number, assumed to be 
conserved. We will not pursue such arguments at this point, preferring to stick to 
the space-time volume form, which is supposed to be available regardless of any 
other considerations. 

5. Reduced elastic constitutive law 

There are different physical and mathematical criteria that can be invoked to 
reduce the general form of the elastic constitutive law (|4.2(l . The first and most 
important of these is the principle of covariance. We will demand that the func- 
tional expressions of the constitutive laws be covariant (that is, tensorial) under 
any change of frame, whether holonomic or not, in spacetime M.. To find the corre- 
sponding reduction in the form of the constitutive law, we commence by noting that 
under a change of frame given by a matrix with entries S l al the frame components 
of the history gradient change according to: 

(5-1) F\ = F%Sl 

where a hat indicates the components in the new frame. Similarly, the frame 
components of the metric vary according to: 

(5.2) 9ij = (S- 1 )US- 1 ) b 3 g ab - 

Suppose at first that ip is a scalar constitutive quantity prescribed by a law on J\f 
of the form 

(5.3) ^ = i>(X A ,g abl FX), 

where F% and g a b are evaluated at X A and K a (X A ), respectively. We suppress the 
dependence on material tensors or constants on Af since these remain inert under 
changes of frame on M. . We demand that in the new frame the value of the function 
xjj remain unchanged, namely 2 : 

(5.4) i> = i>(X A , 9ij ,FX) =4,{X A : {S- 1 T l {S- 1 ))g abl F%Sl). 

This is an identity to be satisfied for all choices of the non-singular matrix {S z a }. 
We may, therefore, choose its entries to coincide with those of F _1 , thus obtaining: 

(5.5) ^> = i>{X A ,F i A Fig ij ), 

where tf> is an arbitrary function. In conclusion, the principle of covariance requires 
that: 

(5.6) V = <MX,C), 

2 A jet bundle-theoretic derivation of 15.41 based on an equivariance property under diffeomor- 
phisms is given in the Appendix. 
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where C = F T gF is the pull-back (by the history k) of the metric g to the body- 
time manifold. 3 

Remark 5.1. An interesting identity: A direct consequence of the reduced 
form of a scalar constitutive law is the existence of an identity satisfied by its 
various partial derivatives. It is not difficult to show, by simply using the chain 
rule of differentiation and effecting a few algebraic operations, that Equation (|5.ti|) 
implies: 

(") F |!ff = 2 (g| g , 




or, for the sake of clarity, in components: 

(5-8) F}°t = 2P-g ia . 

It is quite possible that (perhaps modulo some topological condition) the reasoning 
can be reversed: the satisfaction of this identity would then imply the specific 
dependence of the constitutive law on C, rather than on any other combination of 
F and g. 

If we now repeat the previous reasoning for a space-time tensor-valued function 
(such as the stress-energy-momentum tensor t), the only difference arises from the 
fact that the components of t are affected by the change of frame according to the 
rules for second-order tensors. The final reduced form of the constitutive law (|4.2|) 
is: 

(5.9) t = F [ f (X, C) + A U ® U ] F T , 

where f is a purely body-time tensor which, by analogy, may be called the (deter- 
minate part of the) Kirchhoff stress-energy-momentum. 4 

If the constitutive law happens to be hyperelastic, in accordance with Equation 
(|4.4(l and with the chain rule of differentiation, the determinate part of the Kirchhoff 
stress-energy-momentum tensor is given by: 

dW 

(5.10) f = WC - l +2^. 

It is worth noting that the formal reduction of constitutive laws just effected is 
independent of the fibred nature of the body-time manifold. 

Further restrictions on the constitutive law may arise from a number of causes, 
some of which we briefly discuss at this point, leaving a more detailed discussion for 
later. Since the time-consistency condition can be expressed as C44 = —1, it follows 
that, in any given system of coordinates, there is no need to include a dependence 
on C44. This observation can be used, in combination with other criteria, to further 
reduce the constitutive law. Another restriction may arise from physical reasoning 
whereby, if dissipation is to be ruled out, the vector u should be an eigenvector 



3 tt is interesting to note that this result coincides with the classical one (derived from the 
classical principle objectivity or frame-indifference), by a subtle balancing act between the new 
variables (namely, the presence of a non-trivial metric) and the changes of frame allowed (namely, 
arbitrary changes, rather than just orthonormal). 

Note that this tensor differs from the so-called second Piola- Kirchhoff stress by the absence 
of the determinant of the deformation gradient. In fact, the second Piola-Kirchhoff stress depends 
on the choice of reference configuration and can be more properly understood as a tensor density. 
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of t. Finally, further reductions can be obtained by invoking additional material 
symmetries (such as isotropy) that the material may have. Before returning to 
these matters, we proceed to the formulation of explicit dynamic equations. 

6. The dynamic equations 
A fully fledged dynamical problem calls for the solution of Einstein's equation: 

(6.1) Ein = t, 

where Ein is the Einstein tensor associated with the metric g. The right-hand side 
of this equation incorporates the constitutive law, itself a function of g through its 
dependence on C, according to Equation (|5.9(l . But, since the Einstein tensor is 
divergence-free, we need to impose the following integrability condition: 

(6.2) V-t = 0, 
or, in coordinates: 

(6.3) t a % = 0, 

where a semicolon indicates covariant differentiation with respect to the (torsion- 
free) metric connection. The left-hand side of this equation is ultimately expressible, 
via the constitutive law 1(5.9(1 . in terms of the kinematical variables. Consequently, 
Equations 1(6.1(1 and 1(6.2(1 constitute a total of fourteen partial differential equations 
which, together with the universal constraint 1)3.2(1 . provide the field equations to 
be solved for the fifteen unknowns g a b, K a and A. When the presence of the elastic 
body is likely to affect the gravitational background only slightly, one may use a 
perturbation technique whereby, as a first step, the metric g is kept at a fixed value 
and only the dynamic equations 1(6.3(1 are solved (always in conjunction with the 
constitutive laws and the constraint) for n a and A as functions of X A . This first 
step, often called the fixed-background problem, is usually considered sufficient 
for non-cosmological applications, for which a very slight violation of Einstein's 
equation is certainly tolerable. 

There are many useful manipulations afforded by sundry pull-backs and pro- 
jections of the field equations. The most common one, and one that can readily 
be interpreted physically, consists of resolving Equation 1(6.2(1 into a scalar compo- 
nent on the local four-velocity vector u and the remaining projection on its normal 
three-dimensional hyperplane. Our objective has been to show that at the level of 
generality maintained so far a consistent and complete theory emerges and that the 
kinematical variables, including the Lagrange multiplier, are ultimately obtainable 
as a solution of the field equations supplemented with the equation of constraint 
and the constitutive law. An example of some importance to evidence the physical 
meaning of the Lagrange multiplier is that of dust, that is, a continuous collection 
of material particles without any mutual clastic interaction. In this extreme case, 
we naturally set the constitutive function f in 1(5.9(1 to zero and obtain the following 
residual form of the constitutive law: 

(6.4) t = Au(g)u, 

which can be recognized as the standard stress-energy- momentum tensor of a dust, 
provided one identifies A with the mass-energy density. For the sake of the exercise, 
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if we now enforce the dynamic equation l|6.2(l and take its inner product with u, 
namely, 

(6.5) g(u,V-(Au®u)) =0, 

and if we take account of the constraint i|3.2[l .we readily obtain the scalar equation: 

(6.6) V-(Au)=0, 

which, when subtracted from (|6.5[1 . yields the geodesic equation: 

(6.7) V u u = 0. 

7. A DRASTIC REDUCTION OF THE CONSTITUTIVE LAW 

Given a space-time (scalar-, vector- or tensor-valued) constitutive function ^(n, g, F), 
a constitutive symmetry at the point n £ J\f is an automorphism H of the tangent 
space T n N such that: 

(7.1) $(n,g,F) = $(n,g,FH), 

for all history gradients F at n. The physical meaning of a symmetry is that 
the particular body-time point n is indifferent in its constitutive response to the 
pre-application of a transformation H of its neighbourhood. Notice that this trans- 
formation applies only to the deformation gradient, and not to the various material 
tensors which are the repository of the material properties. It is not difficult to 
show that all the symmetries of a constitutive law at n form a group Ti n , which is 
called the symmetry group of \I/ at n. In the particular case of the constitutive law 
(|5.9(l . a symmetry H must satisfy the following identity: 

(7.2) f(n, C) + AU«U = H f (n, H T CH) H T + A (HU) <g> (HU). 

The first conclusion that imposes itself, by virtue of the independence of the term 
governed by the Lagrange multiplier from the determinate term, is that: 

(7.3) HU = ±U, 

or, in other words, that U must be an eigenvector of H with a unitary eigen- 
value. This conclusion is consistent with the fact that other transformations are 
already ruled out by the time-consistency condition Ij3.2|l . Excluding, moreover, 
transformations involving time reversal, we may say that the eigenvalue is +1. In a 
body-time chart (consistent, as always, with the fibred structure of A/"), a symmetry 
H with components must have the following general matrix expression: 

K 0" 

h T 1 



(7.4) H 



Here, K is a non-singular 3 x 3-matrix and h T is a 3-row. Within the group of 
matrices having this particular form the symmetry group must lie as a subgroup 
and, in principle, it could be as small as the trivial group (consisting of just the 
identity transformation) . 

We now define a particular symmetry group of (|5.9() . which we will denote by 
H e . It consists of all automorphisms H (having, of course, U as an eigenvector 
with unit eigenvalue) that preserve the bundle projection, namely: 

(7.5) tt*(HV) =tt*V, V V e T n J\f. 
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In a coordinate representation, a typical element of this group looks as follows: 
(7.6) H 



I 

h T 1 



with the arbitrary K having been replaced by the unit matrix. Although a perfectly 
valid formulation of elastic materials (according to our definition) can be pursued 
without restricting in any way the symmetry group (beyond the general condition 
(I7.3|0 . most authors seem to implicitly assume that every elastic material must have 
a symmetry group large enough to contain H e . This restriction does not seem to 
arise from any mathematical consideration, but rather from a putatively physical 
reasoning that, by identifying the lack of symmetry of this type with the sensitivity 
of the material to a relative time-shift between neighbouring particles, detects a 
source of energy dissipation. Be that as it may, it will become presently obvious 
that, in order to relate our formulation to such well-established works as those of 
Carter and Quintana or Beig and Schmidt [Q, the condition 



(7.7) 



will have to be adopted. 

We will now investigate the rather severe restrictions imposed by condition (|7.7|) 
on the elastic generic constitutive law ()5.9JI . The term involving the Lagrange 
multiplier needs no further reduction. We are, therefore, left with the identity: 



(7.8) 



f (n, C) = H f (n, H T CH) H 2 



to be satisfied by all non-singular symmetric tensors C (with C44 = — 1) and all 
H G TL e . As a preliminary calculation, we will consider the simpler question of 
finding the restrictions that would apply to a scalar constitutive law at a point 
n G M, namely, we intend to find the most general form of a scalar- valued function 
%jj = ip(C) satisfying the identity: 



(7.9) 



V>(C) = <A(H T CH), 



for all H G TL e . This is certainly a simpler problem than that posed by Equation 
(17.81) . We will carry out the proof of our formula in a local frame and, at the end 
of the process, we will provide an invariant representation of the result. Let, then, 
the matrix expression of C in a local body-time frame (always consistent with the 
fibration) be given by: 



(7.10) 



C 



A 



where condition 13.2|) ') has been used in the form C44 = — 1. The 3x3 matrix A 
is symmetric and positive definite. We now evaluate the matrix product H T CH 
directly from Equations (|7.6|) and 1)7. 1U[) as: 



(7.11) 



H T CH = 



bh T + hb 2 



hh J 



b h 

-1 



Since h is arbitrary, for any given C we may certainly choose: 



(7.12) 



h = b 
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whence it follows that the function tp must depend on its arguments through the 
following peculiar combination: 

"A + bb T " 

-1 



(7.13) 



V>(C)=v 



One may wonder whether further reduction might still arise from other choices of h. 
Nevertheless, it is a straightforward matter to verify that an arbitrary function of 
the form (|7.13|) satisfies the initial identity, so that further reduction is not implied 
by the identity. One may now question the correctness of this formula on the 
grounds that it may not be invariant under changes of (fibre-consistent) body-time 
frames. One way to dispel this fear is to check directly that under a change of frame 
the form of the formula is preserved. A more illuminating alternative is obtained 
by recalling that, given a non-singular symmetric rank 2 covariant tensor Z, we 
can canonically define a rank 2 contravariant counterpart Z , and vice- versa. The 
matrix representations of these entities in mutually dual bases are (as suggested 
by the notation) inverses of each other. It is now a straightforward matter to 
check that the expression A + bb is the matrix representation of the following 
twice-covariant non-singular symmetric tensor at ir(n) £ B: 

(7.14) MCT 1 ))- 1 . 

Note that the inversion inside the projection 7r* is on the four-dimensional manifold 
M., while the outer inversion is on the three-dimensional manifold B. Notice, too, 
that the resulting three-dimensional covariant tensor is symmetric and positive 
definite. 

Having solved the preliminary scalar problem, we are now ready to tackle the 
more delicate tensor identity lj7.8|) . Accordingly, we partition the matrix of compo- 
nents of the body-time symmetric-tensor-valued function f as follows: 



(7.15) 



f(C) = 



M 

T 

n 1 



where the dependence on the point n £ Af has been omitted for simplicity. Since we 
will need to evaluate this matrix function also at H T CH, wc adopt the temporary 
notation of indicating such quantities with a circumflex accent. Thus, for instance, 
the 3x3 matrix M = M(H T CH) is the matrix taking the place of M in the 
representation of f(H T CH). Performing all the matrix operations indicated in 
Equation lj7.8|l . we obtain the following identity: 



(7.16) 



We have at our disposal the degree of freedom of varying h to make this identity 
work. We start by observing that from the upper- left block of the identity we obtain 
that: 

(7.17) M(C) = M(H T CH), 

where we have reverted to the standard notation for function arguments. From our 
previous experience with the scalar law, we conclude that the matrix M (that is, 



"M 


n 




M 


Mh + n 


T 

n 1 


(1_ 




h T M + n T 


h T Mh + h T n + n T h + q 



each of its entries) will necessarily depend on C through the combination A - 
The remaining equations are: 



bb" 



(7.18) 



n(C) = Mh + n, 
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and 

(7.19) g(C) = h T Mh + 2h T n + q, 

where the hatted quantities are to be evaluated at H T CH. The identity l|7.18|l is 
satisfied by the following surprisingly simple form of n: 

(7.20) n = Mb + p, 

where p is an arbitrary function of A + bb T . Indeed, we have: 

(7.21) n = Mb = Mb = M(b - h), 

where we have made use of (I7.17|) and (|7.11|) . Plugging l|7.20l) and (|7.21|) into l|7.18|l . 
we convince ourselves that the identity is indeed satisfied for all h by the function 
n = Mb + p. We should also prove that this is the only possible solution of this 
identity, but we refrain from trying, considering ourselves fortunate to have found 
this solution. As far as the remaining identity, Equation (|7.19|l . plugging into it the 
result H7.20fl just obtained, and remembering that by n we now mean M(b — h), 
we conclude that it is satisfied by adopting the following form for q: 

(7.22) q = b T Mb + 2b T p + r, 

where r is an arbitrary function of A + bb T . It can be verified that, upon coordi- 
nate transformations in the body-time bundle, the quantities M, p and r behave, 
respectively, as a tensor, a vector and a scalar defined in the base manifold. In 
a local frame consisting of U and three vectors perpendicular to it (with respect 
to the pulled-back metric C), we have that b = 0, and the blocks in the matrix 
representing f are precisely M, p, p T and r. 

Notice that the fact that we have restricted the argument C by means of the 
condition C44 = — 1 (a condition that was used repeatedly in the calculations) 
implies that the function f(C) can be determined only up to an additive term of 
the form AU® U, as we already know. This implies an indeterminacy in the choice 
of r. 

Remark 7.1. The condition H7.20[l is automatically satisfied by any hyperelastic 
constitutive law such that W = W(A 4- bb T ). Thus, every hyperelastic material 
whose (scalar) constitutive law is reduced by the action of 7i e gives rise to a stress- 
energy-momentum tensor also satisfying this reduction. The converse, however, is 
not true in general: there exist non-hyperelastic constitutive laws (for the stress- 
energy-momentum tensor) that pass the test. Indeed, there is no a-priori reason for 
M to be derivable from a potential. In fact, if the constitutive law is hyperelastic 
then it turns out that p = and that M is itself derivable (according to the 
hyperelastic prescription |5.1Up ) from a scalar function of A + bb . 

Remark 7.2. We could have refrained from imposing ab initio the condition 
C44 = —1, in which case we would have obtained an apparently determinate com- 
ponent r. This determination, however, would have only been illusory, since the 
general expression for the stress contains a corresponding term affected by a La- 
grange multiplier. The situation is similar to what happens in classical Continuum 
Mechanics when we deal with an internal constraint, such as incompressibility. We 
may use an elastic constitutive law of a compressible material which would provide 
a determinate hydrostatic pressure (via, say, a bulk modulus) for any given pro- 
cess, but this pressure is to be ultimately corrected by the Lagrange multiplier. A 
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different bulk modulus will certainly change the pressure determined by the same 
process, but the Lagrange multiplier will adapt its value to correct the situation. 

Our results so far can be summarized as follows. The stress-energy-momentum 
tensor of a general relativistic elastic material is given by the constitutive law: 

(7.23) t = F [ f (n, C) + AU®U] F T , 

where f is a purely body-time tensor called the Kirchhoff stress-energy-momentum. 
The Lagrange multiplier arises from the requirement that all possible histories pre- 
serve the metric structure of the fibres of the body-time manifold J\f. The strict 
dependence on the pulled-back metric C arises from a requirement of relativistic 
frame indifference. Moreover, although not mathematically necessary, one may re- 
quire on physical grounds that the symmetry group of elastic materials be large 
enough to contain the group Ti. e consisting of all local body-time transformations 
that preserve the body-projection of vectors. In that case, the constitutive law 
is further reduced so that (the determinate part of) the Kirchhoff stress-energy- 
momentum in a local body-time frame Ep, U has the form: 

f = A/ rA E r ®E A 

+ (A/ rA C A4 + /) (E r ®U + U®E r ) 

(7.24) + (M rA C r4 C A 4 + 2C r p r + r) U®U, 
where Af rA , p r and r are functions of the twice-covariant body tensor: 

(7.25) MC- 1 ))- 1 = (C7 rA + C T4 C 4A ) E r <g> E A . 

Finally, we will check that once this drastic reduction has taken place and provided 
that p = 0, the stress-energy- momentum tensor has the 4- velocity vector u as an 
eigenvalue. We check this fact directly by reverting to the matrix notation. We 
start by ascertaining that U is an eigenvector (with respect to the pulled-back 
metric C) of the total Kirchhoff stress-energy-momentum as follows: 



(7.26) 



M Mb 
b T M b T Mb + r + A 



A b 

b T -1 



oi r o 

1/ 1-r + A 



which proves the assertion. The fact that u = FU is an eigenvector of t is now a 
straightforward consequence of the definitions of t and C. In fact, we could have 
convinced ourselves directly of this result by simply adopting in J\f a frame con- 
sisting of U and any three vectors C-orthogonal to it, since under these conditions 
we have b = and, therefore Mb = 0. From these considerations it also follows 
that the (purely space-time) stress-energy-momentum tensor has no mixed (time- 
with space-like) components in a frame consisting of the 4-velocity vector u and 
any three space- like vectors g-orthogonal to it. 

Our purpose in this section has been to show that, except for the presence of 
the Lagrange-multiplicr term, all of the usual assumptions about the form of the 
constitutive law in relativistic elasticity, as pioneered by Carter and Quintana 0, 
are recovered in the present formulation by imposing a particular type of body- 
time symmetry. Specifically, it must be assumed that the symmetry group of the 
constitutive law of all materials under consideration contains the group TC e . It 
appears that, although there is no strictly mathematical reason for adopting such a 
restrictive criterion, there may be relatively strong physical reasons to disregard all 
other materials. From the point of view of the theory of anelasticity that we shall 
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propose, the question is not crucial. The main reason to have delved into the elastic 
realm with such detail has been to make sure that the four-dimensional body-time 
elastic archetype rests on a solid foundation. 

8. Constitutive symmetries: solids and fluids 

We have already introduced and partially exploited the concept of constitutive 
symmetry in Sectional Our purpose in this section is to investigate the presence of 
further constitutive symmetries, assuming that the constitutive law already enjoys 
the standard symmetry introduced in Section^ In other words, we assume that the 
symmetry group TL n of the constitutive law at n € fif contains Ti e as a subgroup, 
as already expressed in Equation (|7.7|l . It is natural, therefore, to introduce the 
quotient group 

(8.1) Q n = n n /n e 

as our object of interest and to call it the reduced symmetry group of the consti- 
tutive law. Naturally, once a chart is chosen in TV, this group is expressible as a 
multiplicative matrix group. A general change of chart (always consistent with the 
projection ir and with the structural group of the body-time bundle) is given by 
four smooth functions Y 1 , Y 2 , Y 3 , y 4 of the form: 

(8.2) Y A = Y A {X A ), A, A = 1,2, 3, 
and 

(8.3) Y 4 = X 4 + y 4 (X A ), A = 1,2,3. 

A moment's reflection reveals that, as far as the reduced symmetry group is con- 
cerned, the function y 4 is irrelevant. In fact, the study of the reduced symmetry 
group in terms of coordinate representations boils down to the study of the sym- 
metries of the matrix function M(A + bb T ). Under coordinate transformations 
in the body-time principal bundle, these symmetries are sensitive only to the part 
embodied in Equation 1)8. 2|l . Another way of expressing these ideas is to say that 
we will investigate the symmetries of the projected constitutive law on the base 
manifold B. For definiteness, but without much loss of generality, we shall concen- 
trate our attention on hyperelastic constitutive laws, characterized, as we know, by 
a single scalar function W of a purely material (three-dimensional) symmetric and 
positive definite twice-covariant tensor Z = A + bb . A (reduced) symmetry of 
such a constitutive law, consists of an automorphism K: 

(8.4) K : T <n) B — » T w(n) B, 
such that: 

(8.5) IU(K T ZK) = W{Z), 

for all symmetric and positive definite twice-covariant tensors Z. 

It is clear that the symmetry groups thus obtained for the same constitutive law 
in two different charts are mutually conjugate, the conjugation being established by 
the gradient of the change of chart. In particular, if the symmetry group relative to 
one chart is unimodular, then perforce it will be unimodular in all charts. Bearing 
this idea in mind we may and shall adopt the standard classification of material 
symmetries of classical Continuum Mechanics. In particular, we have the following 
definitions: 
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Definition 8.1. An elastic material point is a solid if the symmetry group of 
its constitutive law as represented in any chart is conjugate of a subgroup of the 
(Euclidean) orthogonal group. 

Definition 8.2. An elastic material point is isotropic if its symmetry group (in 
any chart) contains a conjugate of the orthogonal group. 

Definition 8.3. An elastic material point is a fluid if its symmetry group is the 
unimodular group. 

Notice that the unimodular group is characterized by the fact that it is the 
largest subgroup of GL(3,R) that preserves all volume forms. In any given chart, 
this implies that the function W depends on its argument (namely, Z) through 
the determinant of Z. In order for this function to be invariant under changes of 
coordinates, therefore, we need to introduce any other volume form (for instance, 
the baryon form) and make W depend on the ratio between the determinant of Z 
and the component of this form. We conclude that the most general constitutive 
equation of a fluid is given by an arbitrary function of this ratio. 

Let the quotient group Q n be known. Then the corresponding original group TL n 
consists of the 4 x 4-matrices of the form: 

K 0" 

h T 1 



(8.6) H 



where K is in Q n and h is arbitrary. In other words, if a constitutive law expressed 
in terms of Z satisfies (I8.5|) . then the corresponding constitutive law in terms of C 
as an independent variable will satisfy a similar equation, namely: 

(8.7) M / (H T CH) = W(C), 

where H is an arbitrary matrix of the form (|8.6(l . The verification of this prop- 
erty and of the fact that these matrices transform in the appropriate way under 
coordinate transformations of the type (|8.2|l and l|8.3[l is a straightforward exercise. 

9. Constitutive isomorphisms and uniformity 

A material body-time complex J\f may enjoy symmetries that go beyond the 
constitutive symmetries of each of its points. These "symmetries" , arising from 
the comparison of the constitutive responses at different points of the body-time 
manifold, confer to it an extra geometrical structure, namely, that of a Lie groupoid. 
In classical Continuum Mechanics it is the time evolution of this entity that allows 
for a rigorous theoretical formulation of the anelastic behaviour characteristic, for 
example, of ideal clasto-plasticity and of bulk growth. In this section we will review 
some of the basic notions of this theory as they apply within the context of General 
Relativity. To emphasize the fact that these notions are quite independent of other 
concepts introduced so far (such as the various reductions of the constitutive law) , 
we will return to the primeval form of a constitutive law for some space-time (scalar-, 
vector- or tensor- valued) constitutive function 'J', namely: 

(9.1) * = *(n,g,F). 

Definition 9.1. Two points rii, € M are constitutively isomorphic if there exists 
an isomorphism P12 between their tangent spaces: 

(9.2) P 12 : T ni N — T n2 M, 
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such that the equation 

(9.3) *(n 2 ,g,F) = #(ni,g,FPia) 

is satisfied identically for all g and for all F in their respective domains of existence. 
Such a P12 is said to be a constitutive isomorphism. For consistency, we will assume 
that the isomorphism P12 respects the unit vector field U, namely: 

(9.4) P 12 U(n x ) = U(n 2 ). 

In local charts consistent with a trivialization of the body-time complex Af, this 
condition is equivalent to P12 having a matrix representation of the form (|7.4|) . 

The preceding definition is the relativistic analog of the classical notion of mate- 
rial isomorphism introduced by Noll f |22U18| 1 It is clear that this definition reverts 
to that of a constitutive symmetry, namely l|7.1|l . if we just identify n\ with n 2 . 
From the physical point of view, we may say that when two points are constitu- 
tively isomorphic they are "made of the same material" . Indeed, when two points 
are constitutively isomorphic, there exist local charts that render the coordinate 
expressions of their constitutive laws identical to each other, and vice-versa. 

Definition 9.2. A material body-time complex Af is said to be constitutively uni- 
form (or, simply, uniform) if its points are pair-wise constitutively isomorphic. If 
the material isomorphisms can be chosen to depend smoothly on both the source 
and the target points, Af is said to be smoothly uniform. 

Due to the assumed fibred nature of the body-time manifold, it is also possible 
to introduce the following notion. 

Definition 9.3. A body-time manifold Af is time-wise uniform if the points within 
each fibre are pair-wise constitutively isomorphic. 

Physically, a time-wise uniform body-time manifold is such that each point in 
the base manifold B preserves its own constitutive nature (its "chemical identity" , 
as it were) as time goes on, without regard as to whether this nature is comparable 
with that of any other point of B. Clearly, 

(9.5) uniformity =>■ time — wise uniformity, 

but not conversely. Because of the "one-sided" bundle nature of the material body- 
time, there is no spatial counterpart to the concept of time- wise uniformity. On the 
other hand, if a body-time complex is known to be time- wise uniform and if for each 
point of the base manifold there exists a local cross-section the points of whose image 
are constitutively isomorphic, then the body-time complex is necessarily uniform. 

Suppose that Af is smoothly uniform. It is clear that if P12 is a constitutive 
isomorphism between the source point n\ and the target point n 2 , then the inverse 
P21 = P12 is a constitutive isomorphism in which the source and the target have 
been exchanged. Moreover, if P12 and P23 are constitutive isomorphisms between 
n\ and tt-2 and between 712 and 713, respectively, then the composition P13 = P23P12 
is a constitutive isomorphism between n\ and 713. Finally, the collection of all 
constitutive automorphisms at each point n constitutes a group (more precisely, 
this is the constitutive symmetry group TL n at that point). It is not difficult to 
prove that the symmetry groups at two different points are conjugate of each other, 
and that the conjugation is achieved by any constitutive isomorphism between these 
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two points. Conversely, given a constitutive isomorphism P12 between n\ and 712, 
the totality V12 of constitutive isomorphisms between these points is given by: 



(9.6) V12 = Pi2H ni = H 



n 2 r 12- 



All these properties taken together confer upon a uniform body-time complex Af 
the structure of a transitive Lie groupoid. 

A local chart (X A ) on Af whose domain contains n and n' naturally induces the 
map O : TAf — > TAf which written as a two-point tensor has the form 



(9.7) 0(n,n') ' } 



dX J 

This gives rise to: 



dX J 



Definition 9.4. A smoothly uniform body-time complex Af is locally homogeneous 
if, for each point n S Af, there exists a chart (containing n) such that the maps it 
naturally induces between tangent spaces are constitutive isomorphisms. 

It is often useful to exploit the transitive character of the material isomorphisms 
to extract a particular point, no say, from Af and to use it as an archetype in 
the sense that the definition of uniformity is equivalent to the following fact: all 
points of Af are constitutively isomorphic to this archetype. Denoting a material 
isomorphism from the archetype no to a point n by P(n), and the corresponding 

set of all such isomorphisms by V(n), we see that choosing a frame {E^} at the 
archetype induces at each point n the collection of frames {Eyi(n) = ■p(n)E^' ) }, 
which is a subset of all the possible frames at that point. This subset is governed 
by the constitutive symmetry of the archetype. 

In the case where the symmetry group is the identity each of the sets V(n) 
contains only P(n). For each such P(n) and frame {E^n)} one may define a 
linear connection V on Af such that (VE J 4)(n) = 0. In this frame the components 
{P A } of 

(9.8) P(n) = P^(n,n )E A (n)®B s (n ) 

are constants, where {B A } is the dual coframe field to {E^} i.e. H a (Eb) = 8g for 
all n. The existence of the V-parallel frame field {E^} implies that the curvature 
of V vanishes. Furthermore, if the {Eyt} are holonomic, i.e. E^ = -g^x in some 
local chart X A , then the torsion of V also vanishes and V is said to be flat. Should 
the torsion of V not vanish, one can say that Af contains a smooth distribution of 
inhomogeneities or dislocations. More generally, the sets V contain more than one 
member and their geometric interpretation is formulated in terms of G-structures 
00IH]- The notion of local homogeneity corresponds exactly to a notion of flatness 
of these G-structures. 

A smoothly uniform body-time complex Af (whether locally homogeneous or not) 
can be described in terms of the archetypal constitutive law ^(g, F) = ^(hq, g, F) 
of a point no m the form: 

(9.9) *(n, g, F) = *(g, FP(n)). 

Assume now that the archetypal constitutive law happens to be hyperelastic and 
that we identify ^ with the elastic energy W. Regarding, accordingly, W as a 
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function of g, F and P(n), we obtain, by virtue of Equation l|9.9|) : 
(9 10) d_K^ ¥ Td_W^T 

[ ' op of ' 

or, invoking (|4.4|l : 

(9.11) — = [~wi + f t { f~ t ] p T . 

which, in components, reads: 

(9-12) S = [-WSf + g ad Fi r b (F- 1 )*} (P" 1 )^. 

The quantity enclosed in square brackets is a purely body-time tensor which, 
following the classical counterpart El E3j j is referred to as the determinate part 
of the relativistic Eshelby tensor . From the above formula, it is clear that its 
physical meaning is related to the elastic energy expended in producing a change 
(or remodeling) of a first-order neighbourhood in the body-time manifold. In view 
of other concepts introduced in this section, this interpretation can be worded as 
follows: the amount of elastic energy required to change the pattern of distribution 
of inhomogeneities. 

If the constitutive equation is given explicitly in terms of an equation such as 
(19.91) . with a specific choice of the field P(n) (with the degree of freedom afforded 
by the symmetry group, of course), it means that we have somehow been able 
to specify the constitutive behaviour at each point of the base manifold for all 
times. In practice, however, starting from an initial space-like Cauchy manifold, 
the constitutive law will evolve at each point in a way that, although always abiding 
by Equation l|9.9|l . is determined by some extra constitutive criterion that allows 
a specific P to be pinned down according to, for example, the local value of the 
Eshelby tensor. Thus, the P-maps function somewhat as internal state variables 
governed by some laws of evolution. 

10. Evolution laws 

Let 1Z denote the canonical right action of R on Af. For each s £ R we have, 
therefore, a fibre-preserving diffeomorphism: 

K s : Af — ► Af 

n i ► TZ s (n) = ns 

(10.1) (X A ,S) h-> {X A ,S + s) 
The induced tangent map: 

(10.2) n„ m : T s U^T ns M 

is, therefore, a non-singular linear map between the corresponding tangent spaces. 
Conversely, given two points, n\,n<i £ Af, such that 7r(nx) = 77(712), there exists a 
unique s £ R such that n-i = n\s and, therefore, a uniquely determined map R12 
between their tangent spaces. Assume now that we are given a time-wise uniform 
body complex Af and that for each pair of points m and ri2 that lie in a given 
fibre 7r _1 (6) the map R42 just defined happens to be a constitutive isomorphism. 
In such a case we are justified in saying that, as far as the point b £ B is concerned, 
the constitutive equation does not evolve in time. If this is the case for all points 
of B, we say that the constitutive equation of Af is non-evolutive. Conversely, if for 
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at least one point b of B and at least one pair of points ni,ri2 G 7r _1 (6) the right- 
action-induced map R12 is not a constitutive isomorphism, we are in the presence 
of material evolution 

As pointed out in Section |2 we can associate with the one-parameter group 
action 1Z the fundamental unit vector field U. Let us denote the Lie derivative 
with respect to U by means of a superposed dot, namely, for a vector field V: 

(10.3) V = LuV=[U,V]. 

Let J\f be a time- wise uniform body-time complex. On each fibre 7r _1 (6) we may, 
therefore, choose a particular point, no(b) say, as a fibre archetype and denote by 
P(n) a smooth choice of constitutive isomorphisms from n Q (b) to the points n lying 
on the corresponding fibre. These linear maps, being two-point tensors, can be 
regarded as vector fields (more precisely, as vector-field-valued covectors at the 
archetype). If the body-time is actually uniform, then a single archetype can be 
chosen. The non-evolution condition can then be stated as: 

(10.4) P = 0. 

Accordingly, an evolution law at a point n € Af will be an expression of the form: 

(10.5) P = *(P,e,...;n), 

where e is the (determinate part of the) relativistic Eshelby tensor (as per Equations 
(|9.11|l or (|9.12(l L and where other arguments could be included. As is the case in 
the non-relativistic counterpart, laws of evolution are subjected to a number of 
formal restrictions, some of which we will presently derive. 

The first restriction on the possible forms of the evolution function stems 
from the uniformity requirement itself. Indeed, if all the points are made of the 
same material, the evolution law of each point c should be obtained as the push 
forward by P(n) of the evolution law at the archetype, namely, a law of the form 
P = <I>o( e o)- As a result of this requirement, we obtain that the evolution law must 
necessarily be of the form: 

(10.6) P = P # (P T eP~ T ). 
Introducing the notation: 

(10.7) Lp^P^P, 

we write the reduced evolution law as: 

(10.8) L Po = * (P T eP- T ). 
The quantity: 

(10.9) L P = PP \ 

can be referred to, by abuse of terminology, as the (general relativistic) inhomo- 
geneity velocity gradient. It is a linear map of the tangent space T n Af into itself. 
The tensor Lp = P _1 LpP is, therefore, the pull-back of the inhomogeneity veloc- 
ity gradient to the archetype. To summarize, the reduction of the general evolution 
law by uniformity arguments leads to an explicit dependence on the uniformity 
map P(n), as given in Equation IjlO.fip . It is worthwhile noting that in a local chart 
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consistent with a trivialization of the body-time complex, the matrix representing 
the tensor Lp has the generic form: 




where L and m are, respectively, an arbitrary square matrix and an arbitrary vector 
of order 3. 

Since the uniformity maps are in general non-unique, and since their lack of 
uniqueness is governed by the constitutive symmetry group, it is to be expected 
that further reductions (due to symmetry) will be possible whenever the symmetry 
group is non trivial. These reductions, to whose derivation we turn presently, are 
of two kinds. 

The first symmetry reduction, called the principle of actual evolution [111 , 
stems from the observation that the condition (|10.4|) is in fact sufficient but not 
necessary for claiming that the constitutive law does not evolve. Indeed, let P(n) 
be a smooth choice of constitutive isomorphisms from no (6) to the points n lying 
on the fibre 7r _1 (6), and let this P(n) satisfy 1)10. 4|l . Consider now a different choice 
of constitutive isomorphisms given by: 

(10.11) Q(n) = P(n)G(n), 

where G(n) is a smooth one-parameter (n) family of material symmetries of the fibre 
archetype such that G(no) = I. Such a choice can be made non-trivially provided 
that the symmetry group is continuous. Clearly, the constitutive isomorphisms 
Q(rt) represent the same material phenomenon as the original P(n), since they 
differ at any point along the fibre just by a material symmetry of the archetype. 
In other words, the choice Q(?i) corresponds to a non-evolving situation along the 
fibre. Nevertheless, Q(n) will in general fail to satisfy condition (|10.4|) . Indeed, 
taking the Lie-derivative of Equation (|10.11|l with respect to U we obtain: 

(10.12) Q- 1 Q = G- 1 G, 

where we have made use of l|10.4[) 5 . It follows then that, as long as the evolution 
function $o gives a result within the Lie algebra of the archetype, there is no 
evolution. The principle of actual evolution, therefore, states that the evolution law 
must have values lying outside of the Lie algebra of the archetype. Two apparently 
different evolution laws whose results differ by an element of this Lie algebra are, 
therefore, to be considered as equivalent. 

Finally, the evolution law must be invariant under the action of the symmetry 
group of the archetype, namely: 

(10.13) * (e ) = G * (G T e G- T ) G 1 + G, 

for all members G of the symmetry group. In this formula, G is an arbitrary 
element of the Lie algebra of the symmetry group of the archetype. The reason for 
the presence of this term should be clear from the principle of actual evolution. 



'Note that G coincides with the ordinary derivative with respect to the (real) fibre parameter. 
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11. Example: an anelastic relativistic fluid 

According to Definition l8.3l and the remarks thereafter, a hyperelastic fluid point 
is completely characterized by a scalar function of the form: 



(11.1) W = W(Vdet(Z)/w), 

where u) is the component of some given 3- volume form #1 on the three-dimensional 
base manifold B: 

(11.2) #1 =uodX 1 AdX 2 AdX 3 

in the chart X T . We recall that, in accordance with the notation of Section |SJ 
Z = A+bb T is the matrix representation of the tensor (7r*(C -1 )) -1 , defined at each 
point of the three-dimensional base manifold B. The above expression \J det(Z) /uj, 
though, is chart independent because it is the ratio between the components of two 
tensor densities of the same type. It is now a straightforward matter to obtain the 
most general expression for the stress-energy-momentum tensor t corresponding to 
a potential W of the form (|11.1|) through the use of Equations (|5.9II and l|5.1U[l . 
The calculation is greatly facilitated by the fact that: 

(11.3) det(Z) = -dct(C), 

as it follows from Equation (|7.1U|) 6 . In other words, for the case of a fluid, it doesn't 
matter whether the space-time volume form is defined from the fully-fledged pull- 
back of the spacetime metric or from its projection on the base manifold. This 
conclusion is also consistent with the fact that the determinant of C automatically 
passes the drastic- reduction test embodied in Equation (|7.9|l . Such considerations 
indicate that the ratio \J det(Z) /u> can be written in a coordinate-free manner. In 
fact 



(11.4) Vdct(Z)/^ = #- 1 [tu(^**l)] 

where tu is the interior derivative on forms, *1 is the volume form of the spacetime 
manifold M. and # _1 is the inverse of the the Hodge map # associated with #1. 
Recalling the formula for the derivative of a determinant, namely, 

(H.5) ^=det(C)C- 

the final form of the stress-energy-momcntum tensor is obtained as: 

(11.6) t = p( v / det(Z)/w) g + Au®u 

as expected, where p is a scalar function of its argument. More precisely, using 
Equation (|4.4J) . we obtain that p is related to W via the formula: 



(ii.7) P =w+ vjggy , 

UJ 

where W indicates the derivative of the function W of Equation 1)11. ljl with respect 
to its argument. 

Assume now that we are given a uniform body-time complex Af modelled after 
a given constitutive law of the type l|11.6l) for the archetype. This means that for 
each point n € fif there exists a linear map P(n) from the tangent space T na M 



^Perhaps the easiest way to convince oneself of this fact is to multiply the matrix 17.101 to the 
whose determinant is clearly 1. 



left by the matrix 



1 
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at the archetype no to T n Af such that the stress-energy-momentum tensor at n is 
given by: 



(11. 



P 



P 



dct(PZP T )/a 



g + Au (g> u, 



(det(P))Vdet(Z)/w g + Au<g>u 



where Z is now evaluated through the pull-back by k of g to n and det(P) is the 
component of the two-point tensor 



(11.9) det(P)(^A 



d d d 
OX 2 A dX 3 A dS 



g> (dX l A dX 2 A dX 3 A dS) 



induced by P. Equivalently, we can write the constitutive equation (|11.8[) in terms 
of C at n as: 



(11.10) 



P 



(det(P))V-det(C)/a 



g + Au £g> u. 



In fact, we can look at this constitutive equation as specifying at each point 
n G M a behaviour of the type (|11.6|l , but with a reference volume two-point form 
defined locally by the formula: 

w 



(11.11) 



dX 1 A dX 2 A dX 3 A dS. 



(det(P)) 

Notice that in this formula u) is the component of #1 at the archetype no, so that 
(|ll.llfl is the contraction of the 4- form #1 A dS (at the archetype no) with the 
two-point tensor 
(11.12) 

-is f d . d . d d 



det(P~ i ) 



\dX l A dX 



dX 3 A dS 



(dX 1 A dX 2 A dX 3 A dS) 



induced by P . 

So as to construct a model of an anelastic fluid fashioned after the archetypal 
constitutive law (|11.1U|) , we will allow the referential volume form to evolve 

according to an evolution law of the type l|10.8fl . For the case of a fluid, the principle 
of actual evolution (introduced in Section ITfj|) . stipulates that the function $o m 
Equation (|10.8(l must not be of the form: 



(11.13) 



K' 

h' T 



where h' is arbitrary and K' is traceless (because the Lie algebra of the unimodular 
group is precisely the algebra of traceless matrices). 

As far as the restrictions placed by Eauation Hll).13fl . we start by noting that the 
determinate part of the relativistic Eshelby tensor is given, according to Equations 
(EU, CDI and EH, by: 



(11.14) 



e = (p - W) I. 



The restrictions just mentioned are satisfied by the following evolution equation 
expressed in matrix form as: 



(11.15) 



p i p 



<f>(p - W)I 
T 
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where <p is an arbitrary scalar-valued function of a scalar argument. Naturally, this 
formula is not form-invariant under (trivialization-consistent) coordinate transfor- 
mations. Nevertheless, the extra terms that would appear affect only the first three 
entries of the fourth row, which is permitted according to the principle of actual 
evolution. In other words, the offending terms would belong to the Lie algebra of 
the group H e . 

At this point, it may prove worthwhile to indicate the non-relativistic counter- 
part of the example we have proposed. Consider a compressible elastic fluid which 
cannot sustain any shear stress, so that its constitutive law consists simply of a hy- 
drostatic stress p determined by the ratio between its present volume V and a given 
reference volume Vq. A possible constitutive law is: p = k(V/V~o — 1), where k is 
positive material constant measuring the elasticity of the material. If, for example, 
this material were to be squeezed into a rigid container of a volume smaller than Vq, 
it would sustain a compression that would persist as long as the material and the 
container are left unaltered. If, however, the reference volume (or resting volume) 
were to begin to decrease with time, the pressure would accordingly decrease. This 
is a volumetric version of the well-known phenomenon of stress relaxation, observed 
to varying degrees in all real materials subjected to a fixed uniaxial extension. The 
rate of decrease of the reference volume may be directly related to the pressure, so 
that as the reference volume approaches the volume of the container, the pressure 
approaches zero and the process tends to fade asymptotically. A possible evolution 
law reads: Vo = hp, where ft. is a positive material constant representing the relax- 
ation properties of the material. In the relativistic picture that we have presented, 
the role of V is played by V '— dot C, while Vo is represented by w/(det P). 
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12. Appendix 

An alternative approach to determine the structure of t is via a variational 
principle. The previous discussion is formulated in terms of the embedding map k 
from the body-time manifold Af to the spacetime manifold (Ai,g). To construct a 
variational principle it proves expedient to employ the inverse / = kT x of k, 

f :M->N 

(12-1) ' _ x 

m I ► TL = K (to). 

Let (£, tt£m, M.) be a bundle with fibre N, and (x a ,X A ) be the coordinates of 
a point in £ with the base point (x a ), i.e. x a — -Kg M (x,X). Thus, / specifies a 
section a of £ with the coordinate representation 



(12-2) 



<j:M->£ 

x a ^(x a ,f A (x)). 



Let J x £ be the first jet bundle of (£, ttsm > -M) with fibre F p at p £ M. and (X A , F A ) 
coordinate points in F p . Take ( / y ab ) to be coordinates on fl p C T p M £g> T p M, the 
space of non-degenerate rank 2 symmetric contravariant tensors at p e M. Let R p 
be the real line at p £ M. and coordinate it by (I). The variational principle below is 
expressed in terms of local sections of a fibre bundle (V, ttx>m , M.) whose fibre over 
p £ M is F p x Sip x R p . Thus, a point in V is coordinated by (x a , X A , T A , j ab , I). 

The inverse metric tensor field G = g _1 on M. and the scalar field AonM are 
used to "prolong" a to the section Jg,\g oi 2? : 



J G ,\o -M ■ V 
(J2 "" ! ' x a ~ (x a ,f A (x),d a f A (x),G ab (x),X(x)) 

where G ab = G(dx a ,dx b ) = g~ 1 (dx a , dx b ) are the components of g -1 with respect 
(x a ) on M. 
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The action functional has the form 

(12.4) S[T] = [ pu> 

Jv[M} 

where T = Jg.\g and u) is the naturally induced volume 4-form 

(12.5) w = ; 1 =dx 1 A dx 2 A dx 3 A dx 4 

y/- det[ 7 Qb ] 

on V and p is a 0-form on T>. 

We consider here media where the 0-form p is constructed out of functions on T> 
including the components of a 1-form fi = pA{X)dX A on N C T>: 

(12.6) p = W(X,F,i) - \l (p A (X)p B (X)T£T* 1 ab + l) 

where the scalar field W on T> characterizes the mechanical properties of the 
medium. It will be shown how the 1-form fi, map / and the metric tensor g 
are used to recover the 4- velocity vector field uonM. 

12.1. Equivariance under diffeomorphisms. Not all scalar fields on T> are 
admissible candidates for W because not all choices lead to spacetime covariant 
variational equations on A4. The dependence of p on r y ab and T A must be such 
that (J7g,ac)*(p w ) is a spacetime coframe independent 4-form on M.. The object 
{Jg,\&)*u is such a 4-form on M. by definition so only {Jg,\&)*P need be consid- 
ered. Geometrically, the local covariance of (Jg,\&)*P on M. may be formulated in 
terms of an equivariance property of Jg,\g acting on p under the group of diffeo- 
morphisms on M. since these can be used to induce local coframe transformations. 

Let ip be any local diffeomorphism on M. with tp = tp^ 1 . Admissible choices for 
p are defined to satisfy 

(12-7) <P*(JG,x<rTp = (J&,x<r)*P 

where 

(12.8) Jl x a = J^a^x^a. 

Let (x a ) be the coordinates oi m E A4 and let p be the 0-form 

(12.9) p(x a ) ee p (x a , f A (x),d a f A (x), G ab {x), X(x)) 
on M. and note that the local coordinate expression of l|12.7(l is 
(12.10) 

p (<p a (x)) = p (<p a (x), (<p*f A )(x), d a (p*f A )(x), (^G) ab (x) 7 (ip*X)(x)) 
= p (<p a (x), (f A o <p)( x ), d a (f A o <p){x), {G cd d c rd d ^ b o <p)(x), (A o p)(x)) 

Equation (|12.7|l may be used to derive an identity analogous to (|5.4|) at mo € M. 
by a restriction to diffeomorphisms for which mo is a fixed point. At m = mo, 
where <p a {x ) = x%, using (|12.9|) and l|12.10|) it follows that 

(12.11) p (xg, /o 4 , (d a f A )o, G a », Ao) - P (x a J A , (S-^Udtf^o, S* S b dn G c d , A ) 

where Sg, = (d b ^ o <p)(xo), = f A (x ), (d a f A ) o = d a f A (x ), G a » = G ab (x ) 
and lo — 1(xq)- Hence, by choosing <p appropriately 1)12. ll[l holds at all points in M 
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for any S% and, following the same arguments given earlier, may be used to reduce 

the dependence of p: 

(12.12) 

p {x a J A (x) 1 d a f A (x),G ab (x),X(x)) = p (x a , f A (x), G ah d a f A d b f B (x),\{x)) 
for some function p. 

An identity analogous to (|5.8|) may be obtained directly from l|12.7fl by using a 
1-parameter family of local diffeomorphisms <p e generated by the vector field Y on 
U C M: 



(12.13) 
(12.14) 



</?o(w) — m V to € W C M., 
d 



Using iflTTjl it can be seen that 
d 



(12.15) 



de 



H>\ (Ja,\<r) P 



s=0 
d 

de 



dx a ' 



6 = 



where ip £ = (p £ 1 and so (|12.1U|) leads to 

.. dp 
dx a 



£ Y p = y^ + £ Y /^ + a o (£ Y /^ dp 



(12.16) 



(£ Y G 



dX A 



dT A 



where each partial derivative of p is implicitly evaluated on Jg,\& & n d Y a 
dx a (Y). Furthermore, (|12.9|) yields 



^p=y a &+^f A ^+^(d a f M dp 



(12.17) 



fcx d P 



dx a ' ~ 1J dX A 
+ £ Y (G afc ) 

which subtracted from l|12.1(j|) leads to 



dT A 



Q^ab 



r \ dp 



(12.18) 



d a Y b d b f 



dp 



2G a 



dp 



dT A dj cb 
The previous equation holds for any Y a and so 



(12.19) 



d b f A j^ = 2G ac 9(3 



q t a dl cb 
and furthermore, using i|12.6fl it follows that 



(12.20) 

on the section Jc\cr of T>. 



. BW dW 



dT A 



dy 



cb 



12.2. The stress-energy tensor and normalization constraint. The action 
functional (|12.4|) is varied with respect to G via variations of the section Jg,ac 
to obtain the stress-energy tensor T and is extremelized with respect to / and A 
via J7g,ac to obtain the equation of motion for / and the normalization constraint 
g(u,u) = -1. 
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(12.21) 



Let r e = Jg s ,\p be the 1-parameter family of sections given by 

x a i * (x a , f A (x) : d a f A (x), Gf(x), X(x)) 

where G E is a 1-parameter family of contravariant metric tensors on M. and Go = 
G. The covariant stress-energy tensor T on M, dual to the contravariant stress- 
energy tensor t, of the medium is defined by 



(12.22) 



de 



S\T E 



le=0 



SG ab T ab * 1 



M 



where * is the Hodge map associated with the volume 4- form *1 = TqU! on A4, 
T a b = T(d a ,dh), uj is given in (|12.5(l and 



(12.23) 

with 

(12.24) 



de 



= 6G 



ab 



d 



6=0 



8G 



ab 



dG 



ab 



de 



r= 



ab 



is the variational vector field associated with JG et \cr. Using 
(12.25) ^(det[ 7 ah ]) = det[ 7 afc ] 7cd <5 7 cd , 

where 7 ab 7 b c = 5° (the Kronecker delta), leads to 

dp 



(12.26) 



T„, = T* -2 



Q-jCib 



Plab 



Although T has a similar structure to the stress-energy tensor employed in recall 
that here / has rank 4 and TqW involves both space and time derivatives of /. 
Now let T £ = JG,A e cr be the 1-parameter family of sections given by 

Jg,X s (7 :M^V 

(12 27) 

{ ■ ' x a ^(x a J A (x),d a f A (x),G ab (x),X £ (x)) 

where X £ is a 1-parameter family of scalar fields on M. and Ao = A. Demanding 
that 

d 



(12.28) 

for all A variations 

(12.29) 

with 

(12.30) 

leads to 
(12.31) 



de 
de 



S[T e 



6 = 



= sx 







dl 



sx = 



dX e 
de 



6 = 



HA(f)p B (f)df A A*df B = -*l 

where d is the exterior derivative. The vector field u is constructed from the 1-form 
f*fi = [iA{f)df A and the metric g as 

(12.32) a(u) = G(/>,a) 
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where a is any 1-form on M. and, using equation H12.31|) . u satisfies the normal- 
ization condition 



where u a are the coordinate components of the f-form dual to the vector field u. 
The equations of motion for / can be obtained in a similar fashion, or derived by 
setting the divergence of T to zero. 

We refer to (|12.34|) as the stress-energy tensor for a simple relativistic hyper- 
elastic medium. In this article the words elastic and hyper- elastic are used purely 
by analogy with concepts in classical (non-relativistic) Continuum Mechanics. Since 
k% involves time derivatives this is clearly an abuse of terminology. 

The derivation of the form of the hyper-elastic stress-energy tensor l|12.34|) based 
on (jl2.7|l is often referred to as compatibility with the principle of local covariance 
on A4. The subtlety here is that (J7g,ac)*(pw) must be a 4— form on M. when 
constructed out of tensors on the body-time complex M and the metric tensor 
g on M. This covariance is the general spacetime analogue of classical material 
objectivity (for an account of the classical theory see |22p. 

In general the stress-energy tensor may not arise from an action of the form (|12.4(l 
but still may be written in the form (|4.2(1 for some symmetric degree 2 tensor f and 
where n = f(x) = n^ 1 {x). 
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